Coarsening and accelerated equilibration in mass-conserving heterogeneous nucleation 
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We propose a model of mass-conserving heterogeneous nucleation to describe the dynamics of 
ligand-receptor binding in closed cellular compartments. When the ligand dissociation rate is small, 
competition among receptors for free ligands gives rise to two very different long-time ligand-receptor 
cluster size distributions. Cluster sizes first plateau to a long-lived, initial-condition dependent, 
"metastable" distribution, and coarsen only much later to a qualitatively different equilibrium one. 
Surprisingly, we also find parameters for which a very special subset of clusters have equal metastable 
and equilibrium sizes, appearing to equilibrate much faster than the rest. Our results provide a 
quantitative framework for ligand binding kinetics and suggest a mechanism by which different 
clusters can approach their equilibrium sizes in unexpected ways. 

PACS numbers: 82.60.Nli,02.30.Hq,05.70.Ln 



I. INTRODUCTION 



The binding of multiple particles to specific nucleation 
sites is a key process in many physical and chemical set- 
tings. The formation of droplets, condensates on aerosols 
[l|, |2] , and crystals Q is often triggered by the pres- 
ence of impurities or boundaries, in a process known 
as heterogeneous nucleation [3, Q. Heterogeneous nu- 
cleation also occurs in cell biology during the assembly 
of sickle hemoglobin @, /3-amyloid fibers 0, Arp2/3 
complex-mediated actin nucleation Q , and probably dur- 
ing clathrin-coat assembly Q . Within biochemical appli- 
cations, ligand-receptor binding can also be viewed as a 
particular paradigm of heterogeneous nucleation, where 
multiple ligands bind to a single receptor akin to an im- 
purity seed in solid-state nucleation. 

Viewed through this lens, nucleation is ubiquitous in 
cell biology. Indeed, receptor loading levels control a va- 
riety of biochemical reactions, from viral entry to cell sig- 
nalling. The chemical stoichiometrics involved in ligand- 
binding events however, may limit the maximum num- 
ber of ligands a receptor can hold to about a dozen. 
For example, hemoglobin can bind at most four oxygen 
molecules [lO| , virus-cell fusion occurs after a small num- 
ber of cell surface receptors bind to a viral protein [ll| , 
and cell-signaling is initiated after a certain number of 
phosphates bind to specific enzymes p^ . This is in con- 
trast to most physical and chemical systems where ag- 
gregation of an unlimited number of particles can lead to 
the emergence of macroscopic structures. 

An even more critical feature of nucleation in cellu- 
lar settings is the small system volumes involved and, 
as a consequence, the presence of a finite number of 
monomeric ligands driving the nucleation process. In 
the small volumes encountered in cells, ligand produc- 
tion and degradation are often slower processes than at- 
tachment and detachment to receptors, allowing certain 
ligands to be depleted 13| . Because there is no source to 
replenish the free ligand concentration, receptors in con- 
fined, isolated systems compete amongst themselves for 
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FIG. 1: A heterogeneous nucleation process in which ligand 
monomers bind only to seeds. Here, Ns = 6 seeds (open 
hexagons) are available to bind M — 30 initial monomers 
(filled dots). 



the finite pool of free monomeric ligands, as depicted in 
Fig. [TJ For simplicity, throughout the remainder of this 
paper, the terms "monomers" and "ligands" will be used 
interchangeably, as will "seeds" and "receptors." 

The dynamics of mass-conserving homogeneous nucle- 
ation has been well-studied in the context of Becker- 
Doring equations [l^ - [l6| . In this work we study its het- 
erogeneous counterpart, relevant for ligand-receptor ki- 
netics in biology. While heterogeneous nucleation has 
been well studied, most theoretical treatments focus on 
computing equilibrium partition functions for nucleation 
with specific forms for the free energy of monomer asso- 
ciation (T7l - [l9| . Many other approaches focus on either 
the molecular details and geometry of an individual clus- 
ter particle [20| . or on the asymptotic dynamics of even 
more coarse-grained continuum size distributions (2lj | . In 
many ap plic ations, a constant source of monomers is also 
imposed [22| . Here, we will instead consider the dynamics 
of a system with a total fixed number M of monomeric 
ligands (bound and unbound) and a fixed number Ns 
of receptor seeds. Each receptor can bind at the most 

monomeric ligands according to the spatially-uniform 
mass-action equations we describe in the next section. 

Two qualitatively different cases are analyzed. In Sec- 
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tion III, we first consider irreversible binding, where the 
detachment rate is strictly zero so that once attached, 
monomers cannot detach from clusters. Irreversibility 
leads to a loss of ergodicity since only a fraction of the 
possible cluster configurations will be sampled during the 
dynamics, while many others will never be visited. As 
a result, the final "quenched" or "metastable" cluster 
size distribution is not an equilibrium one and depends 
strongly on initial conditions. Ergodicity is restored in 
the case of a non-zero detachment rate, where all pos- 
sible cluster configurations are eventually sampled and 
where the cluster sizes approach an equilibrium distribu- 
tion, independent of the initial configuration. Reversible 
binding, in the limit of small unbinding rates, is in ana- 
lyzed Section IV. 

Several studies of homogeneous nucleation have shown 
the existence of long lived metastable states followed 
by final equilibration, or "coarsening" , to a very differ- 
ent cluster size distribution. These results were found 
for nucleating systems driven by an infinite supply of 
monomers and while allowing clusters to grow without 
bound [IH • O^'^ results for mass-conserving het- 
erogeneous nucleation show a similar coarsening behav- 
ior. The steady-state cluster distributions arising from 
the irreversible and reversible dynamics are very dif- 
ferent from each other, especially in the limit of small 
particle numbers and even in the case of vanishingly 
small detachment rates. In the latter case, when un- 
binding is very slow compared to binding, relaxation to 
the true equilibrium cluster size distribution occurs over 
the long time scales associated with unbinding. As a re- 
sult, cluster concentrations reach long-lived metastable 
plateaus that depend on initial conditions and that can 
be closely approximated by results obtained from con- 
sidering irreversible dynamics, as treated in Section III. 
Only at longer times, after monomers start unbinding 
in appreciable numbers, does this metastable size dis- 
tribution "coarsen" and cross over to the true equilib- 
rium one. While the metastable and equilibrium clus- 
ter distributions are generally very different, we find a 
surprising result: for certain sets of parameters, special 
cluster sizes have identical metastable and equilibrium 
concentrations. For these clusters, the equilibration pro- 
cess appears to be dramatically accelerated. In Section 
V we find the exact mathematical relationship leading to 
the apparent fast coarsening where certain clusters reach 
equilibrium concentrations well before the rest. Our re- 
sults are a consequence of total mass conservation, and 
do not arise in the case of receptors binding an unlimited 
supply of free ligand monomers. Finally, in the Conclu- 
sions, we discuss implications and future extensions of 
our work. 



II. MASS-ACTION EQUATIONS 

To begin our analysis, we consider a model of hetero- 
geneous nucleation for M well-mixed monomeric ligands 



binding sequentially [2^ to any of the iVg uniformly dis- 
persed ligand seeds, neglecting fragmentation and aggre- 
gation that do not involve monomers, since they have 
been treated in other contexts [l^ [2^ . We also assume 
that each seed can accommodate at most N monomers 
due to stoichiometry constraints and consider the mean- 
field mass-action equations for the number of clusters 
Cfe(t) of size k, where < k < N. Here, fc = indicates 
"naked seeds" , with no bound monomers, and k — N 
saturated ones, where no further binding is possible. In 
general, monomer attachment and detachment rates from 
a cluster of size k can be explicitly fc-dependent. Specific 
forms for pk and have been used to describe cooper- 
ativity and the nucleation of clusters of various shapes 
and in different dimensions. For instance, pk ^ k^l"^ and 
constant arc typically used to model 2D nucleation of 
circular droplets when monomer binding is not diffusion- 
limited [26j . In our work, we assume that while detach- 
ment is independent of the number of particles in the free 
monomer pool, the attachment process depends only on 
how many monomers remain unbound. The "Becker- 
Doring" equations for Cfe(t) can thus be written as 

Co = -po'm{t)co + qici, 

Ck = -Pkm{t)ck ~ QkCk + Pk~im{t)ck-i + qk+iCk+i, 
Cn = -qNCN + PN-im{t)cN~i- 

(1) 

Here, pkm(t) and qk represent monomer attachment and 
detachment rates, respectively. The attachment rate is 
proportional to m{t), the number of free monomers avail- 
able for binding 

N 

m{t) = M-J2kck{t). (2) 

fc=i 

Note, that although these equations are written assuming 
finite particle numbers, they can also describe concentra- 
tions, given a normalizing reference concentration. The 
quantities M, N, and Ck{t) therefore need not be inte- 
gers. We assume a typical initial condition where all the 
mass is in the form of monomers, m(t = 0) = M,co(t = 
0) = Ns,Ck>o{t = 0) = 0. The other constraint particles 
must obey is that the total number of seeds must be Ng 
at all times, regardless of cluster population levels. We 
thus impose 

N 

iV. = ^c,(t), (3) 

which is satisfied by the system in Eq. [TJ using the given 
initial conditions. For clarity, and because we will be 
referring to these equations often, we rewrite our mass- 
action equations for the simplified case of uniform at- 
tachment and detachment rates pk = p and qk = q- This 
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approximation might be most relevant for modeling nu- 
cleation and growth of linear filaments where there are 
always only one or two ends on which monomers can bind 
or detach from. Rescaling time by the attachment rate 
p, we find 



Co = —m{t)cn + eci, 

Ck = -m{t)ck - eck + m{t)ck-i + ec^+i, (4) 

cn = —ecN + m(t)cN-i, 

where e = q/p. In the context of Eqs. |4l irreversible 
binding corresponds to e = and reversible binding to 
e > 0. In the following, we shall be interested in the 
difference in behavior between e = and e — > 0+ . As we 
shall see, the presence of a vanishingly small detachment 
rate £ — > 0+ can lead to qualitatively different cluster size 
distributions compared to those obtained in the purely 
irreversible case at e = 0. 



dcp 
dr 

dcfc 
dr 

dcAf 

dr 



-Poco, 



(6) 



= Pn~iCn-i- 



Our goal is to find the rescaled time r* corresponding to 
the rescaled time at which monomers are irreversibly de- 
pleted: TO(r*) = M — ^Cfe(T») = 0. The quenched, 
steady-state cluster size distribution is thus found by 
evaluating the cluster concentrations at t*: CA;(r*) = cj!. 
Eqs. [5] are linear and can be solved by using Laplace 



transforms. Upon defining Ck{s) - 
find that Cfc(s) satisfy 

SCO - Ns = -poCo, 



e ''''ck{T)dT we 



III. IRREVERSIBLE BINDING 

We first consider the strictly irreversible binding limit 
in the general framework of Eqs. [1] where there is no de- 
tachment and qk = i). Two possibilities arise. For AI > 
NNs, there is an excess of available monomers. Given 
the irreversible nature of the dynamics, all Ns seeds will 
be fully occupied by N ligands, leaving M — NNs free 
ones. In this case, we expect steady state solutions 
to yield €^{1 — )■ co) = Ng, Ck^jLNit — >■ oo) — and 
m{t oo) ~ M — NNg. We shall call this the excess 
monomer limit. In the other case of M < NNg there 
are not enough monomers to fill seeds to capacity and a 
non-trivial steady state will arise. Here, we expect the 
existence of a finite time at which the pool of free 
monomers is depleted, so that m(t*) = 0. At this time, 
the final cluster distribution is the one frozen at t^, , since 
no further attachments nor any detachments are possible. 
We shall call this the excess seed limit. Since the quantity 
M — NNs will play an important role in our analysis, we 
introduce the monomer excess parameter a = M/NNg, 
so that values of ct > 1 correspond to the excess monomer 
case, while a < 1 describes the case of excess seeds. 

To determine the final cluster size distributions in both 
cases, first note that Eqs. [T] (or Eqs.|3]) are nonlinear due 
to the constraint on m(t) involving Ck(t) via Eq. O If 
qk = however, all terms on the right hand side multiply 
m{t). We can make analytic progress by dividing by m(t) 
and defining a rescaled time r according to 



SCk = -PkCk +Pk^lCk~l, 

scn = pn~iCn-i, 
which yield the solutions 



(7) 



■prfe-l 

ULois+Pj)' 



cn{s) 



■prAf-l 



sUf=o\s+Pj) 



.(8) 



To simplify our analysis, we restrict ourselves to uniform 
intrinsic attachment rates Pk = P and use units of time 
such that p ~ I. The dynamics arc now described by 
Eq.dlwith e = 0. The solutions represented by Eqs. [5] 
thus simplify to 



I 



cn{s) 



N, (s + 1)^+1' A^, s(.s + l)^' 

which can be inverse Laplace transformed to yield 



(9) 



Ck<N{T) 

Ns 



fc! 



Ns 



N-l 



1-E 



(10) 



These results obey the constraint in Eq.O The value 
of can now be found by using Eqs. [TU] in the mass 
constraint Eq. [2] and imposing the condition ?7i(r*) = 0. 
After some algebra, this condition yields 



dr 
dt 



N 



i{t) = M -^kck{t). 



(5) 



k=l 



Eqs. [T] can now be written as 



rfe-"' 
NT{N) 



+ {N-T,) 



r(A^,T,) 
NT{N) 



1-cr. 



(11) 



As mentioned above, we expect a finite solution only 
in the excess seed [a < 1) case. For excess monomers 
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(cr > 1) we do not expect a finite time at which monomers 
are depleted. Indeed, the left hand side of the above ex- 
pression is positive and monotonically decreasing, imply- 
ing that Eq.[TT]will have a finite, real solution only in the 
excess seed case, for a < 1. When the initial monomer 
number M is increased, and a decreases past unity, the 
root T, diverges since all binding sites on the seeds are 
eventually occupied and further depletion of monomers 
can never occur. The quenched concentrations in this 
case are described by m{T^ — >■ oo) = M — NN^ and 
Cfe(T* — >■ oo) = NgSk^N, indicating that all seeds are filled 
to capacity for ct > 1, as expected. 

As a nontrivial example of the seed excess case, cr < 1, 
we can numerically solve Eq.[Tl] for N — 10, M — 30, 
iVs = 8 and a = 3/8 to obtain t, = 3.75248 and 



Co 




Ns 




C2 






= 2!^ 


C4 






= ¥^" 


C6 




Ns 




C8 




Ns 


^8!^ 


ClO 






= 10!^ 



0.088031, 
= 0.206596, 
= 0.145455, 
= 0.0487661, 
= 0.0095373, 

(12) 

It can be explicitly verified that these solutions obey 
Efe=oCfc = Ns. Finally, in the case of large iV, asymp- 
totic analysis of Eq. [TT] gives 



0.02346, 

■ = 0.165168, 

■ = 0.193812, 

■ = 0.0909699, 

■ = 0.0228743, 
^* = 0.0035788. 



Cl 






Ns 


= r*e 




C3 


r! 


-T 


Ns 


= 3!^ 




C.5 


r! 


-T 


Ns 


= 5!^ 






tJ 


-T 


Ns 


= 7!^ 




cg 


r! 


-T 


Ns 


^9!^ 





r* 



M e "-A's 
Ns mV2^ 



eM 



N 



(13) 



which allows us to find analytic approximations to the 
final quenched values = Cfc(r*) in this limit. The full 
irreversible dynamics are illustrated in Fig. [21 In order 
to find approximations for the maximal concentrations 
of clusters Ck{t) of size k < N, we note that from Eg.fTUl 
they occur at the rescaled time t^^x = k. To find the 
corresponding real time tmax we insert Eqs. [lU] into the 
scaling relationship Eq.[5] so that 



= M - NsJ2 




dr 
dt 



Equation 1141 can be numerically integrated to find t{T). 
However, for large N the first sum on the RHS can be 
approximated by its — >■ oo limit 



N-l 1. _ 

r'-e 



k=l 




FIG. 2: (a) Numerical solution to Eqs. [4] with M = 50, 
= 6, A's = 5, and e = 0. Since monomers are in ex- 
cess (cr = 5/3 > 1), all clusters except cjv vanish at long 
times. This plot is indistinguishable from the one plotted us- 
ing e = 0.0001, and is qualitatively similar to what would be 
found for a constant free monomer concentration m{t) = M. 
(b) The numerically computed (colored dots) maximal cluster 
concentrations Cfe(tmax) = fc*e~'"/fc! and corresponding times 
imax. The approximation tmax ~ N^'^ In [M/{M — Nsk)] and 
the corresponding Cfc(fniax) are also shown by the dashed and 
solid curves, respectively. 



In this same limit, the second sum is 0(1/A). Ea.[T4lcan 
therefore be accurately approximated by 



dr 

— w M - NsT, 
dt 

which can be analytically integrated to give 



(16) 



Oil/N). 



(15) 



i(r) « In [Af/(il/ - A,r)] . (17) 

We plot W = t(w) = m « A^-i In [M/(M - A^.fc)] 
and the associated Cfc(iinax) = k'^e^'^/kl in Fig.[2][b) as a 
function of k. The approximation for tmax is shown by 
the dashed curve in Fig. [IJb) and is extremely accurate, 
especially for small k where maxima are reached before 
appreciable accumulation of larger clusters invalidate the 
approximation t M — NsT. 

In the next section we will analyze the nucleation pro- 
cess when successive monomer detachment is allowed. 
The question will arise as to how closely the irreversible 
nucleation results found here are followed in the case of 
a vanishingly small, but non-zero, detachment rate e. As 
we shall see, our reversible results will closely mirror the 
irreversible ones in the limit e — > 0"*", only in the excess 
monomer case, when seeds are saturated with ligands. 
In the excess seed case on the other hand, dramatic dif- 
ferences between reversible and irreversible binding arise, 
even as e — 0"*". Only very special parameter choices will 
lead to the rare matching of reversible and irreversible 
dynamics for specific clusters. 



IV. REVERSIBLE BINDING 

In this section we find the equilibrium cluster size dis- 
tributions allowing for positive detachment rates qk > 0- 
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We start by finding the equilibrium cluster concentra- 
tions cj^ = Ck(t ^ oo) by setting dck/dt = in Eqs. [T] 
Due to reversibility, initial conditions are irrelevant. Af- 



ter defining m®'' = M — X^feLi ^'^'k^^ '^^ '^T' 



be written as a function of Cq'^ and m"'^ as follows 



n, = l 13 



(18) 



This expression can be used in the mass constraint of 
Eq.[2]and the total cluster number constraint in Eq.[3]to 
find two equations for the two unknowns Cq'^ and m"^ 



w 2 ■ 



N=40 




1/2 

G=M/iN,N) 



k=l 



fe=i 11^=1 % 



iCqifc 



(19) 
(20) 



These equations can be solved by substituting the ex- 
pression for Cq'^ in Eq.[20]into Eq. ll9l and determining m°'i 
numerically. Again, computations are greatly simplified 
by restricting our analysis to uniform attachment and de- 
tachment rates pk = p and qk = q, respectively. Further 
nondimensionalizing time in units oip~^ and introducing 
e = q/p, Eqs. [TSl become 



„cq _ cq 



(21) 



cq 



where z 

in Eq.[3] yields c, 



The fixed seed number constraint 
= Nsiz ~ l)/(z^+i - 1) so that by 
substituting c^z^ into Eq. [2]wc find an equation for z: 



£Z 

iVJV 



<j]{z-l){z 



N+l 



l) + z 



N+2 



(22) 



-(l-M/iV).^+i + ^ = 0. 



Eq. [22] determines the numerical value for the normal- 
ized cluster fugacity z. In the small detachment limit 
e — > 0+, once more, the two limits of excess monomers 
and excess seeds naturally arise. In the excess monomer 
case, cr > 1, the real root of Eq.[22]can be found as the 
inner solution of a singular perturbation (27| where the 
largest power of z multiplies e so that z = £~^{M — 
N^N) + Ns{M - NsNy^ + 0(e). Inserting this ap- 
proximation for z into the seed constraint Eq.|3] wc find 
Co"' w [e/{M - N^N)]'^ + 0{e), which yields 



-N-k 



{N,N)N-^ {a - 1) 



N-k 



0{e 



N-k+l 



). (23) 



Thus, equilibrium concentrations all vanish as 
0{e^~'^) except that of the maximum cluster k ~ N 



FIG. 3: Values of the normalized cluster fugacity z deter- 
mined from the real root of Eq. 1221 Fugacities for N = 
3,10,40 are plotted in the limit e/{NsN) = 10"^. A good 
approximation to z is given by Eq. [21] and is shown by the 
dashed curves. The region a = M/{NsN) > 1 corresponds to 
the case of excess monomers where ~ £^~''', as shown in 
Eq. [23 



which asymptotes to w iVg — 0{e). This qualitative 
behavior is expected in the excess monomer case when 
nearly all available binding sites arc occupied and only 
nearly fully occupied seeds survive. In particular, when 



e 



+ oq 



-^k^N 



and 



TVs. This result is identi- 



cal to what was found for the strictly irreversible case of 
e = in the previous section. In the excess monomer case 
(cr > 1) thus, all clusters will be filled to capacity in the 
case of vanishingly small detachment rates, regardless of 
when the limit e 0"*" is taken. 

In the opposite case of excess seeds, cr < 1, monomers 
are depleted before all binding sites on all A^s seeds can be 
filled, leading to finite concentrations c^''. Interestingly, 
the excess seed limit further separates into two sub-cases. 
From our numerical analysis of Eq.[22]we find that z > 1 
for 1/2 < cr < 1, implying c^j^-^ > c^'^ and larger cluster 
sizes tend to be favored. On the contrary, for cr < 1/2 
we find z < 1 so that < c^. In this case, there are 
too few monomers M for larger clusters to persist and 
smaller cluster sizes are more populated. For a range of 
values of cr near 1/2 we find that the approximation 



2 - 



1 - 



24 



N + 2 



cr 

2 



-,1/2 



(24) 



and the associated c?. 



are highly accurate. 
= 1/2 the monomer 



cl^iz)z'' 

Note that at the special point cr 
fugacity z = 1 and all equilibrium concentrations 
Ns/{N + 1) are equal. The behavior of the root z of 
Eq. [21]as a function of the monomer excess is plotted for 
cr < 1 in Fig. [3] The analytic approximation (Eq. [24| is 
also indicated by the dashed curves. 

Our analysis thus far does not provide insight into 
how the equilibrium state is reached. As discussed ear- 
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lier, when e 0+, we expect binding to occur in 
a nearly irreversible manner over intermediate times, 
yielding metastable cluster size distributions. Repeated 
monomer detachment and reattachment become signifi- 
cant only after much longer times, of the order ~ £~^, 
allowing redistribution of mass into equilibrium clusters. 

To find the metastable cluster size distribution, we 
make the ansatz that [t) can be approximated by set- 
ting the detachment rate e = at intermediate times. 
We may thus neglect detachment and use the results ob- 
tained for irreversible binding up to ~ £~^, beyond 
which detachment effects may become appreciable, both 
in the excess monomer and excess seed cases. 

Fig. shows the full time dependence of Cfc(t) in 
the reversible, excess monomer case where a = 5/3 > 
1. Here, as expected, both c|.'^^ and Ck<N(j oo) = 
Cfc<Ar vanish as e — > 0+. In this case, the dynamics is 
not appreciably affected by the onset of detachment and 
there are no dramatic behavioral crossovers originating 
across the time scale. Reversible and irreversible 
dynamics thus coincide at all time scales in the e — )• 0+ 
limit if (T > 1. In particular, the initial rise in co<:k<N{t) 
is determined by the monomer loading process and is 
independent of the detachment rate e. Setting e = 
and using our results from the previous section we can 
numerically compute the time scale imax over which the 
cluster size distributions peak. We have verified that as 
e — )■ 0^, the full dynamics arising from the reversible 
binding process in not appreciably different from that 
obtained in the irreversible binding limit. In summary, 
when monomers are in excess, cr > 1, the differences 
between c^'' and vanish in the e — >■ 0+ limit, ergodicity- 
breaking is not apparent and seeds are always filled to 
capacity. 

We now consider reversible binding in the cr < 1 case, 
where there are more receptor seeds than initial free 
monomers. As in the ansatz made in the excess monomer 
case, we assume that at least up to time scales of order 
£~^, the dynamics can be approximated as an irreversible 
binding process, where e = 0. By following the full dy- 
namics in Eqs. |4] we verify that at intermediate times, 
the metastable concentrations Ck(t) approach levels ap- 
proximated by the final ones Cfe(T*) = = NsT^e~'^' /k\ 
reached in the irreversible case when e = 0. Following 
the evolution of Eqs. S] beyond timescales £~^, we find 
that these metastable concentrations eventually coarsen 
towards a qualitatively different, equilibrium distribution 
defined by ^'^(e ^ 0+). 

This qualitative difference in cluster size distributions 
is noticeable only in the excess seed limit when cr < 1 
and illustrates ergodicity-breaking at e = 0, as clearly 
shown in Fig. g] In Fig. |l{a), where tr = 3/8 < 1/2, 
seeds are in such strong excess that monomers are quickly 
depleted and c^,^,^ < c^. Fig.SJJb) plots c^^ as a function 
of e, found from numerically solving Eq. [22l Note that 
the values of c'f^is 0+) differ from the intermediate 
ones approximated by the frozen distribution c^. The 
latter are indicated by the colored dots. Figures |3Jc) and 




FIG. 4: Ergodicity breaking occurs only when seeds are in 
excess, (a) Numerical solution to Eqs. [4] with e — 0.0001, 
M = 30, iV = 4, and = 20 (a = 3/8). In this strong 
excess seed case, cr < 1/2 and both cl > c^^i and c^"* > 



(b) c^i^ as a function of e. Note that even as e 



0^ 



are different from the metastable values (colored dots) 
found from setting e = in Eqs. [D (c) Cluster concentrations 
Ck{t) for £ = 0.0001, M = 30, iV = 4, and TVs = 10 (cr = 
3/4). In this weak excess seed case, a > 1/2, both < d^^i 
and c^f^ < c^-^. (d) Again, ergodicity-breaking arises since 
c'j!^{e — s- 0) T^- cX{e — >■ 0). In all plots, the free monomer 
concentration m{t) and the number of naked seeds Co{t) can 
be reconstructed from the constraint conditions and are not 
explicitly shown. 



(d) are the analogous plots but for 1/2 < cr < 1, where 
'^l+i '^l- When seeds are in excess, the crossover to 
equilibrium is clearly observable over the coarsening time 
scale tc ~ £~^. 



V. APPARENT ACCELERATED 
EQUILIBRATION OF SPECIFIC CLUSTERS 

The general qualitative behavior described in the pre- 
vious section is that when seeds are in excess (cr < 1), the 
full heterogeneous nucleation problem exhibits dynamics 
occuring over two time scales. The first is oi t ~ 0(1) 
and corresponds to monomer attachment rates, while the 
second coarsening time scale, tc ~ £~^, is associated with 
the monomer detachment rate. In general, cj^ ^ c^. 

However, upon fine tuning relevant parameters, we find 
special values of cr < 1 and N where up to two specific 
cluster sizes k can have nearly equal values of quenched 
and equilibrium concentrations {cj^ ~ c^) in the e — > 0+ 
limit. These clusters quickly reach their equilibrium con- 
centrations on a short time scale independent of e. Math- 
ematically, the sizes k that are subjected to this rapid, 
apparent equilibration can be found by setting = c^f^: 
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7V+1 



1)' 



(25) 



respcc- 



where and z are determined by Eqs. llll andr 
tivcly. 

Fi gure [STa) shows the relative difference (c^ 
for fixed TV = 6 and cr as a function of the discrete cluster 
size k. Generally, we find that many values of < (t < 1 
give rise to at least one value of k at which quenched clus- 
ter concentrations equal equilibrium concentrations. In 
the above example, N = 6 and a = 0.35633, and clusters 
of size fc = 4 (red arrow) quickly quench to their equi- 
librium values. Fig. [Sjb) plots the numerical solution to 
Eq. Hfor e = 10^1", iV = 6 and ct = 0.35633. For sim- 
plicity, we have plotted only C3(i),C4(i), and c^it). Note 
that even though c| = (dashed line), C4{t) does suffer 
a small transient perturbation due to the rearrangement 
of all other clusters Cfc^4 at time t ^ temporarily dis- 
turbing the balance of C4(i). Figs.[5l[c) and (d) illustrate 
the behavior for a — 0.86293. Here, Fig. [SJc) predicts 
that ci quickly reaches its equilibrium value. Fig. [5{d) 
explicitly plots co(i),ci(t), and C2{t) for a — 0.86293. 
Figs. \5la.) and (c), also suggest that c| « over a 
wide range of values of a. For the finite processes we 
have considered, we find that at most two sizes k, out 
of TV, can exhibit accelerated equilibration, provided a 
and N are precisely tuned. We expect the qualitative 
aspects of these results to hold when binding and / or un- 
binding rates are weakly cluster size-dependent, allowing 
our analysis to apply in scenarios of weakly cooperative 
ligand- receptor binding 24| . 



VI. CONCLUSIONS 

In this work we have analyzed a simple and mathe- 
matically tractable model of heterogeneous nucleation to 
describe ligand-receptor binding in closed systems, such 
as within cells or organelles. A complete analysis in 
terms of the parameters {M,N,Ns) in the e — >■ 0+ limit 
shows how dramatically differently the system behaves 
in a monomer-rich environment compared to a seed-rich 
one. In the latter case, when binding sites outnumber 
initial monomers (cr < 1), we find that after an initial 
transient ^ tmax, cluster densities first approach c^., ap- 
proximating the quenched concentrations when e = and 
all free monomers have been depleted. This long-lived 
metastable distribution eventually coarsens to a very dif- 
ferent equilibrium distribution c^,'^ at much later times 
tc ^ 1/e. Surprisingly, when parameters (cr and N) are 
finely tuned, it is also possible that « c^'^, for par- 
ticular fc— clusters, resulting in much shorter coarsening 
times. We find that clusters of up to two specific sizes 
may appear to reach their equilibrium concentrations by 



t (1). Our results have general implications for ligand- 
receptor kinetics and suggest practical ways of tuning 
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FIG. 5: (a) The relative difference (c^ — c°^)/c'j^ as a func- 
tion of discrete values of k for a = 0.08, 0.35633, 0.92. (b) 
Selected cluster concentrations for a = M/{NsN) = 0.35633, 
where c^^4 ~ c^'l^. Note the small transient in C4,{t) near 
t ^ . (c) (cl — cj^)/c^ plotted as a function of k for 
a = 0.08,0.86293,0.92 These values indicate that for a = 
0.86293, the concentration ci{t) quickly reaches its equilib- 
rium value, (d) The corresponding concentration plot show- 
ing just co{t),ci{t), and C2(t). These plots also indicate that 
for a — 0.08, the concentration C2{t) experiences accelerated 
equilibration. 



(M, N, Ns) in experiments to accelerate the equilibration 
of specific clusters by stabilizing their metastable sizes. 

A number of extensions of our analysis can be further 
investigated. Certain forms for cluster size-dependent 
attachment and detachment rates, pk and qk, can be in- 
corporated into the analysis. For example, if pk ''^ fc, 
certain products and sums in Eqs. [THl and [^D] can 
be analytically expressed or approximated to derive vari- 
ations to Eq. [21] and the associated concentrations cj^. 
Furthermore, for small numbers of clusters, the mean- 
field results derived from the Becker-Doring equations 
may deviate from the expected cluster size distributions 
arising from fully stochastic simulations [13] ■ We expect 
our mean-field results to be qualitatively valid when clus- 
ter correlations are included in the dynamics. A careful 
quantitative investigation of stochastic effects will be in- 
cluded in future work. 
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